Spectrum of the Frobenius-Perron 
operator for systems with 
stochastic perturbation 

Andrzej Ostruszka^ and Karol Zyczkowski^ 

^Instytut Fizyki im. Mariana Smoluchowskiego, Uniwersytet Jagiellonski, 
ul. Reymonta 4, 30-059 Krakow, Poland 

^Centrum Fizyki Teoretycznej PAN, 

al. Lotnikow 32/44, 02-668 Warszawa, Poland 

e-mail addresses: ostruszk@if.uj .edu.pl, karol@cft.edu.pl 

May 17, 2001 



1 Introduction 

Consider a discrete dynamical system / acting on a phase space Q. Time 
evolution of a classical density p (a probability distribution function on Q) 
is governed by the Frobenius-Perron operator M |jl|, 

p' {y) = {M p){y) = I 5{f{x) — y)p{x)dx, where / p{x)dx = 1 . (1) 
Jn Jn 

Its spectral properties are of considerable interest, since they influence prop- 
erties of the correlation function IQ, which may be measured in a physical 
system. Spectrum of the FP-operator was calculated for classical dynami- 
cal systems including tent map, Bernoulli shift, baker map ^ ^ and the 
so-called "four legs" map 0. However, results obtained may depend on the 
choice of the function space FP-operator operates in. Such a case was explic- 
itly demonstrated by Antoniou et al. [§, in which two different spectra for 
a certain random dynamical system were found. One might ask therefore, 
which of those constructions lead to physically meaningful results. 

In all physical situations the system under consideration is inevitably 
subjected to a stochastic perturbation which we usually want to reduce as 
much as we can. The idea we want to present here is that the presence 
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of a perturbation of a small amplitude may choose one particular space of 
functions and thus define a "physical" decomposition of the FP-operator, 
without ambiguities related to the choice of the space of the eigenstates. Note 
that the natural invariant measure (so-called SRB measure 0) of a classical 
dynamical system may be defined in a slightly analogous way, as the unique 
invariant measure stable with respect to stochastic perturbation. Moreover, 
it is instructive to recall the quantum mechanical problem of calculating 
corrections to degenerated energy levels via perturbation theory. A priori 
none of the basis in subspace of degenerated states is distinguished. However, 
one may distinguish a certain basis by applying an arbitrary perturbation and 
later decreasing its amplitude to zero. 



2 Model system 

We analyze a model system introduced in |jlO| which allows one for an exact 



representation of the FP-operator describing dynamical system with noise as 
a finite dimensional matrix. The construction proceeds as follows. 

For simplicity we consider a discrete dynamical system / acting on the 
interval VL = [0, 1) and subjected to an additive noise (with periodic boundary 
conditions) 

Xn+l = f{Xn) + in (mod 1) , (2) 

where ^i, ^2? • • • are independent random variables fulfilling 

stationarity "^(Cn) = '^(0 5 
zero mean (^„) = , 
finite variance {inim) = cr'^^mn ■ 

We choose for a noise for which the probability of transition from x to y 
{Vix,y)=PiO) is 

homogeneous Vi^x, y) = V{x — y), (3a) 
periodic V{x,y) = V{x + l,y) = V{x,y + 1), (3b) 

N 

decomposable V{x,y) = ''^^ AirUr{x)vi{y), (3c) 

l,r=0 

where A = {Air)i^r=o,...,N is a real matrix of finite size {N + 1) of expansion 
coefficients and {ur)r=o,...,N, ('yi)i=o,...,Ar are two sets of linearly independent 
real valued functions. 
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Figure 1: Transition kernel VN{f{x),y) for the map defined by with e = 
and the noise strength characterized by = 20, which gives a ~ 0.0694. 



All these conditions are satisfied by the trigonometric noise 

VN{x,y) = C^cos''{n[x-y]), (4) 

where N is even and the normalization constant — 



assures 



r(Af/2+l/2) 

lo '^N{,x)dx = 1. The noise strength may be parametrized either by or by 
its variance IITOtl 



N 



27r2 



^'(A^/2 + 1) 



27r2 ^ 

k=N/2+l 



2^ h2 ] 12 27r2 ^ P' 

k=l 



(5) 



where stands for the derivative of the digamma function |jTT|. Asymptot- 
ically cr 



AT — 1/A^ — > 0. In this case Air = (^)Sir, while ui(x) = vi(x) 



sin'(7ra;) cos^^'(7rx). 

For a deterministic system, the action of the FP-operator is described by 
(|l|). In the presence of a stochastic perturbation this equation needs to be 
modified. 



p'iy) = / 'P{f{x),y)p{x)dx. 



(6) 



Due to the decomposition property (|3|) we may write 



N 1 

p'iy) = ^^ir Ur{f{x))vi{y)p{x)dx 
^ r rl 



r=0 



Ur{f{x))p{x)dx 



My) 



(7) 



where Vr = '^f=oAirVi. Thus any initial density is projected by the FP- 
operator into the vector space spanned by the functions Vr- Eventually, we 
can represent the action of this operator by a matrix of size + 1 acting on 
the vector of expansion coefficients. Writing p(x) = J2i c^i'^ii^) we obtain 



p(y) = ^MyWk = ^My) d 



kl O-l 



kl 



where D^^ = J dx Uk{f {x)) vi{x) . It is worth emphasizing that the matrix 
representation of the FP-operator for the system with noise is finite dimen- 
sional, while the FP-operator connected with the deterministic case acts in 
an infinite dimensional space. 



3 Examples 



In [|r2| it was conjectured that for continuous dynamics the spectrum of the 
system with noise should tend to the spectrum of the operator describing the 
deterministic system, if the amplitude of noise tends to zero. This correspon- 
dence is based on solving eigenvalue problem for the Foker-Planck equation 
- which a priori is not a simpler task than calculating the spectrum of the 
Frobenius-Perron operator. In our approach we can find approximation of 
the exact spectrum just by matrix diagonalization. The larger dimension of 
the matrix D, the smaller variance cr^ of the noise, and the better the ap- 
proximation of the spectrum of the deterministic system. In order to show, 
how the spectrum of the Frobenius-Perron operator could be approximated 
by our procedure we calculate it for a simple one-dimensional map. 

First of all we have to consider stability of the spectrum if the system is 
subjected to a stochastic perturbation. To address this problem it is useful to 



introduce a notion of essential spectrum ||13| , It is a part of the spectrum 
contained in the smallest disc of radius r such that all eigenvalues outside 
it are isolated and of finite multiplicity^. Blank and Keller have shown |[I3[ 



"'^For linear piecewise expanding ID maps it was shown in |T5[ that r = hm„^oo VTA\f" 
where A/ is the smallest absolute value of the derivative of /. 
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Figure 2: Map given by Eq. (|) for e = 1/20 
together with its Markov partition denoted by dotted hnes. 

that for piecewise expanding maps (without periodic points, for which the 
derivative /' is not defined) subjected to any local perturbation, eigenvalues 
outside essential spectrum are close to those of the deterministic system. 

We analyze a family of piecewise expanding maps which have isolated 
eigenvalues outside the essential spectrum 



and subject it to the perturbation (^) (non local in sense of fl^). Since 
the absolute value of the slope is constant and equal to 2, the radius 

of the essential spectrum equals 1/2 independently of the parameter e. For 
e = there are two eigenvalues with modulus 1 (±1) - the system separates 
into two subsystems Qi = (0, |) and ^2 = (|, 1) which exchange positions 
with each other under every iteration of the map. This fact is related to the 
presence of A2 = — 1 in the spectrum. If we increase e then the modulus of 
the negative eigenvalue will decrease. However, for sufficiently small e the 
subleading (with second largest modulus) eigenvalue is still located outside 
the essential spectrum so we expect to approximate it with our procedure. 
For certain values of e we may construct Markov partitions of Q, write down 
the corresponding stochastic transition matrix T and find analytically its 
subleading eigenvalue A2 which determines the rate of convergence to the 
equilibrium. Fig. g presents map (^ for e = 1/20 together with its Markov 




2(x + e) + 1/2 X < 1/4 -e 
3/2 - 2(x + e) X E [1/4 - e 
2(x + e) - 3/2 x> 3/4 - e 




(9) 
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Figure 3: Subleading eigenvalue as a function of noise width a for map defined 
by plotted for e = 1/20 and e = 1/84 together with the deterministic limit 
represented by horizontal lines. 

partition. In this case 

/ 1/4 3/4 \ 

1/2 1/8 3/8 

1 0' ^ ' 
\ 1/3 2/3 / 

and Aa = -(1 + y5)/4 ^ -0.809, while for e = 1/84 we have As ~ -0.9638. 
Fig. 1^ presents the second eigenvalue A2 of the FP-operator for the system 
with noise obtained by numerical diagonalization of (A^ + 1) dimensional 
matrix D as a function of the noise width a for two different values of e. With 
decreasing noise strength a the second eigenvalue tends to its deterministic 
counterpart (represented by dashed horizontal lines). For e = 1/20 a fair 
approximation is obtained already for N = 200 {a ~ 0.02), whereas for 
e = 1/84 we need smaller noise to obtain satisfactory results. Fig. ^ shows the 
motion of four eigenvalues with the largest moduli in the complex plane as the 
noise strength is decreased. Convergence of the eigenvalue outside essential 
spectrum to its deterministic counterpart is much faster then these of the two 
other eigenvalues belonging to the essential spectrum (one may even question, 
whether they at all converge to the deterministic value). In the same way 
we expect that eigenvectors tend to their deterministic counterparts. Fig. | 
presents the second eigenvector obtained from the matrix representation of 
the FP-operator in the noisy system (e = 1/84) for a ~ 0.02 (dashed line) 
and from the transition matrix of the Markov partition. 

Eventually, we can compare the values of the subleading eigenvalue of the 
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Figure 4: Plot of four largest eigenvalues in the complex plane as a function 
of the noise strength for e = 1/84. The darker color, the smaller noise, which 
changes from cri ^ 0.1 to (T2 ~ 0.02 (A''i = 10, A^2 = 200). Crosses represent 
the eigenvalues of the transition matrix T with the largest moduli, while 
dashed circle of radius 1/2 represents the essential spectrum. 
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Figure 5: Eigenvector corresponding to the second eigenvalue obtained from 
the matrix representation of the FP-operator of noisy system (e = 1/84) for 
a ~ 0.02 (dashed line) and from the transition matrix T of the deterministic 
system (sohd line). 
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spectrum with the decay rate of the autocorrelation function. The correlation 
function of any function h{x), defined by a time average 



1 ^ 

Ch{n) = hm - y h{x{t + n))h{x{t)) (11) 
t=l 

is called autocorrelation function in the simplest case h{x) = x. For an 
ergodic system it can be expressed by an average over the phase space with 
respect to the invariant measure p{x) 

C{n) = I dx p{x) x{n) X , where x(n) = /"(x) . (12) 
Jo. 

In the case considered, in which the dynamics is composed of two parts — 
deterministic and stochastic (see (0)) — one has to average additionally over 
different realizations of the stochastic perturbation. Let us consider first the 
one-step correlation function 

C{1) = l{x{f{x)+0)p{x)dx, (13) 

where the angle brackets denote averaging over different realizations of the 
noise (^). Employing property (^ we can write 



C(l) = J V{f{x),x')xx'p{x)dxdx' = ''^^ J Uk{f{x))vk{x')xx' p{x) dxdx 



Vk{x')x'dx'^ (^j 



Uk{f{x))x p{x) dx] = V ■ u , (14) 



where we have defined v and Vl as vectors of the integrals, J dx'vk{x')h{x') 
and J dxp{x)uk{f {x))h{x) for = 0, . . . , A^. In an analogous way we obtain 
the correlation function of any observable h for longer delay times 



Ch{n) = y 'P(/(xi),X2) . . .P(/(x„),x„+i)/i(xi)/i(x„+i)p(xi) 

= vD"~^u, n>0, (15) 

where (A^ + l)-dimensional matrix D is already defined in (H). Thus one 
may expect the correlation function to be composed of a + 1 exponentially 
decaying modes determined by the moduli of the eigenvalues of the matrix 
representation of the FP-operator. An example of such a situation is pre- 
sented in Fig. ^. The autocorrelation function is calculated by choosing 10^ 
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Figure 6: Correlation function for map with e = 1/84 and width of noise 
a ~ 0.045. The inset presents C'{n) = \C{n) — C(oo)| in semilog scale. 
Exponential fit gives logC" ~ — 0.1198n while the subleading eigenvalue of 
the matrix D gives log IA2I = —0.1201. 



points according to the uniform distribution in Q and evolving them by sys- 
tem / given by (|^) with e = 1/84 (and the noise characterized by = 50). 
The value of the decay exponents obtained from the best exponential fit 
(see inset in Fig. ^ agrees up to 2% with the one we get from the second 
eigenvalue of the matrix representation of the FP-operator. 

Note that the observed decay of correlation was studied for an ensemble of 
points distributed uniformly in Q, that is according to the invariant measure 
of the system. This seems to be more natural approach than in the case 



recently analyzed by Weber et al. |I6| , |T^, in which the initial density was 



determined by the selected eigenvector, while the observed decay rate was 
governed by the corresponding eigenvalue. 



4 Conclusions 

The aim of this paper is twofold. We point out that by introducing a stochas- 
tic noise into a deterministic system and later tending with its strength to 
zero on may distinguish a "physically" important part of the spectrum of 
the associated FP-operator. In this way we suggest to define a "physical" 
spectrum of the classical map, as this robust with respect to stochastic pertur- 
bations. Moreover, we provide a method for approximation of the spectrum 
and eigenvectors of the FP-operator by applying a suitable noise decompos- 
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able in the sense of (|3^). This technique seems to be more justified from the 
physical point of view than just truncating the infinite dimensional matrix 
representation to some finite dimension e.g. (thus effectively introducing 
some perturbation), since one constructs a system with noise of known prop- 
erties and can decrease the amplitude of the perturbation in a controlled way. 
Another methods of introducing noise into deterministic systems to analyze 
their spectral properties were recently used in |l^ . 



We demonstrated that the eigenvalues of the FP-operator located out- 
side the essential spectrum are robust not only against local perturbations 
as proved in [O, but also against non-local perturbations of the form (M. 



Moreover, these stable eigenvalues have a direct physical meaning: they de- 
termine the rate of the exponential decay of correlation in the system. Thus 
our approach of analyzing dynamical system with a stochastic perturbation 
of a variable strength allows one to identify the physically important part of 
the FP spectrum without mathematical ambiguities of selecting a space, in 
which this operator acts. 

On the other hand it would be interesting to analyze spectral proper- 
ties of dynamical systems in presence of a stochastic perturbation fulfilling 
properties ( |3^ ) but different than (^) studied in this paper. We expect the 
eigenvalues of the FP-operator not belonging to its essential spectrum to be 
weakly dependent on the specific form of the probability distribution P(^), 
but this conjecture requires further verification. 

We would like to thank W. Slomczyhski and P. Pakohski for fruitful 
discussions. A.O. gratefully acknowledges financial support from Subsy- 
dium FNP 1/99. This work was supported by a Polish KBN grant no 
2 P03B 072 19. 
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